R=120;
L=100E-6;
C=33E-12;
tau=L/R;
dt=tau/100;
t=0:dt:6*tau;
u1(1)=0;
u1(2)=0;
u2(1)=0;
u2(2)=0;
lt=length(t)
u1=randn(1,lt);
for count=2:1:length(t)

u2(count)=2*u2(count-1)-u2(count-2)+u1(count)-2*u1(count-1)+u1(count-2)-(tau/dt)*(u1(count)+u1(count)*((dt*dt*)/L*C));    
end;
plot(t,u2/3)